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Re-studying the non-Fermi liquid one-particle Green functions (NFLGF) we have extended the 
work of A. Balatsky (Phil. Mag. Lett. 68, 251 (1993)) and L. Yin and S. Chakravarty (Int. J. Mod. 
Phys. B 10, 805 (1996)), among others. We use the moment approach of W. Nolting (Z. Phys. 255, 
25 (1972)) to compute the unknown parameters of the NFLGF's in the framework of the Hubbard 
model. The zeroth order moment requires that our one-particle Green functions describe fermionic 
degrees of freedom. In order to satisfy the first order sum rule a renormalization, 7 7^ 1, of the 
free electron mass is called for. The second order sum rule or moment imposes a relation between 
the non-Fermi liquid parameter, a, the Coulomb interaction, U, and the frequency cutoff, ui c . We 
have calculated the effect of the mass renormalization factor, 7, on some physical quantities, like: 
1) the correlated momentum distribution function, n c (k), close to the effective chemical potential, 
at T = 0; 2) the superconducting critical temperature, T c ; and 3) the superconducting critical 
interaction, A cr , and compare them with analytical results found in the literature. Also, we have 
calculated, for the first time, the isotope effect, a', for non-Fermi liquid systems, which reduces to 
a' — 1/2 (the BCS result) when a — > 0. As a new case of non-Fermi liquid systems, in Appendix 
|X| we have studied two inequivalent coupled Hubbard layers (ICHL) for which we calculate the 
one-particle spectral functions on the layers and perpendicular to them. We discuss the new features 
which appear due to the shift in the two effective chemical potentials and propose some experiments 
to detect the features found from our expressions. 

Pacs numbers: 74.20.-Fg, 74.10.-z, 74.60.-w, 74.72.-h 



I. INTRODUCTION 

The unusual normal state properties in high temperature superconductors (HTSC) Q| have lead to the idea that 
the usual Fermi liquid theory is non longer valid. As a consequence of this, several phenomenological models g|| have 
emerged with the purpose of explaining the strange metallic behaviour of the normal state phase and the behaviour 
of T c vs doping in the HTSC. We should like to point out that the origins of a non-Fermi liquid ground state for a 
strongly correlated material (U > 2D, where U is the local repulsive Coulomb interaction and 2D is the width of the 
free band) in dimensions higher than one, namely, d > 1, is an issue that has not been solved analytically up to the 
present moment. The main arguments which have been cited as responsible for the failure of Landau theory are: 

1. The high critical superconducting temperature is attributed to the CuOi planes, but is well known that in one 
or two dimensional systems superconductivity (or any spontaneous symmetry breaking) is suppressed due to 
fluctuation effects for T^0; 
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2. 



The superconducting phase is very close to the long range magnetically ordered phase, and the exchange inter- 
action J in this strange metal can be too strong; 



3. The concentration of holes x, or the carrier number, is too low; 

4. The Coulomb interaction may be too strong so the adiabatic assumption in Landau's theory may no apply; 

5. For optimally doped materials, the electrical resistivity on the plane, p a b cx T down to T c . For a Fermi liquid, 
Pab °c T 2 . There is not consensus regarding the origin of the anomalous in-plane transport. Two broad classes 
of theories attribute the anomalous behaviour either to singular forward or large-momentum scattering [Q. 
Furthermore, Boltzmann transport equation gives an expression for the magnitude of the resistivity in terms of 
band parameters and a mean-free path between quasi-particle collisions. At low temperatures this expression 
suggests a mean-free path which is much larger that the lattice constant, as in conventional metals. However, 
at higher temperatures the resistivity smoothly increases to large values, suggesting a mean-free path which is 
smaller than the lattice constant, implying the breakdown of a quasi-particle picture 0. 

6. The overdoped materials exhibit a ratio of resistivities, p c / p a t ^ f{T) as in standard anisotropic materials. On 
the other hand, underdoped materials exhibit a divergent out-of-plane p c as T decreases, even if p a b is metallic. 
The c-axis data are not as universal between different cuprate families. 

7. The Hall coefficient, Rh(T), in HTSC in the normal state shows the striking NFL behaviour ||. Rh{T) follows 
a Curie-Weiss type T-dependence and \Rh(T)\ » l/\n e\ for T — > in the underdoped compounds, where 
n is the carrier number. Moreover, Rh(T) > for hole-doped compounds and Rh{T) < for electron-doped 
ones, although each of them has similar hole-like Fermi surface. 

8. The NMR relaxation rate, 1/Ti cx T° for these materials. Remember that for a Fermi liquid, 1/Ti cx T. 

9. In conventional metals, one observes a Drude peak in the optical conductivity at u) = 0, which broadens but 
persists for high temperatures. In contrast, in strongly correlated metals most of the spectral weight is in broad 
features at high frequencies. Furthermore, the Drude peak only exists at low temperatures. 

10. In conventional metals, the thermopower is linear in temperature, has values less than fcs/e w &7pV/K and has 
the same sign as the charge carriers. In strongly correlated metals it can have a non-monotonic temperature 
dependence, can change sign and have values of the order of ks/e. 

11. The isotope exponent, a\ in HTSC is unconventional in different respects. Optimally doped samples show a 
very small a' of the order of 0.05 or even smaller, in contrast to the BCS value of 1/2. This unusually small 
value in connection with the high value of T c lead to early suggestions that the pairing interaction in the HTSC 
cuprates might be predominantly electronic in origin with a possible small phononic contribution. However, 
to complicate things a little bit more, the isotope exponent factor, a', also shows an unusually strong doping 
dependence, reaching values of 1/2, in some cases higher, in the underdoped, T c reduced, compounds. 

Among one of the models to explain the anomalous properties of the HTSC we mention the one of Varma et al 
M, which tries to fit the linear resistivity data of the HTSC by means of the quasi-particle lifetime of the form 
r w l/[T(£fc — fj,)], where T is the absolute temperature and p is the chemical potential. This model is strictly 
phenomenological. There is a second model introduced by Anderson || , whose starting hypothesis is the occurrence 
in two dimensional (2 — d) systems of a state similar to the one from the one dimensional (1 — d) Luttinger liquid [Q. 
The point of view adopted by us is that the HTSC are materials which can be treated by the Luttinger liquid model, 
specially in the underdoped regime. Being so, we explore some consequences of this assumption. 

In such a non-Fermi system the one-particle spectral function A(ej:,uj) = —1/tt Im[G{e^ 1 u) + i0 )], where G(k,uj) 
is the one particle Green function, satisfies the scaling property 

A(A ££ ,A^) = A"- 1 x A(e £ ,u>) , (1) 

where a, the non-Fermi liquid parameter, is given by the condition < a < 1/2. In the case of a Fermi liquid 
system a — 0. We mention that Wen ||, as we do in the paper (Eq. (|l5|)), has shown that the exponent a is not 
universal, since it depends on the coupling constant between the electrons. In the case of the Hubbard model, a 
depends on the local Coulomb interaction. The model that we are going to study by means of the scaling relation 
(Eq. (0)) has poles and branch cuts. This scaling behaviour is a generalisation to dimensions d > 1 of the fermionic 
propagator from one dimension systems (d = 1). It has been showed by different authors that the one-particle 

Green function (1PGF) can be expressed as 



G ^= {»-e jm^e AW- . -«c<a,<« (2) 

v (T,fc'' V p,k' 

where E a rp\ £ represents the spin (charge) energy spectrum of the carriers, u) c is frequency cutoff introduced to maintain 

the dimension of the 1PGF correct and g(a), e 1 ^ are renormalization factors introduced in order to recover the usual 
properties of the Green's functions. Outside the interval [— lu c ; lj c ] the Green function will have the usual 1 /u) behaviour 
as to — > oo. 

Eq. (||) was studied in some detail by Yin and Chakravarty Let us mention that Eq. (Q) was used by 

Chakravarty and Anderson jljj to obtain an interlayer tunneling Hamiltonian (ILT) . Just recently, Chakravarty, Kee 



and Abrahams |15| ] have used the interlayer tunneling Hamiltonian to explain the experiments of Basov et al |16| on 
the c-axis optical sum rule in some HTSC. The structure of the paper is as follows. 

In Section [n] we show that a band renormalization parameter 7 is necessary to satisfy the first order sum rule of 
Nolting [fl7| . The presence of 7 seems to us very natural since when correlations are present the free electronic band 
has to be renormalized too. In the same section, we find that a, w c and U are closely related. The parameter 7 is going 



to play an important role in the following development. In Section III we have recalculated the following quantities 



1) the momentum distribution function, n c {£j:) close to the effective chemical potential, namely, ~ Me//! 2) the 



superconducting critical temperature, T c ; and 3) the superconducting critical interaction strength, A cr . In Section IV 
we present the calculation of the isotope exponent or coefficient, a' , using the non-Fermi liquid one-particle Green 
function (NFLGF). In Section we present our conclusions and the outlook of our line of work. 

While the semi phenomenological theory of Anderson || suggests anomalous exponents, i.e., a^O, a satisfactory 
derivation of them and other details await further theoretical development ]T^ | . In Appendix |A] we have found a novel 
application of our theoretical treatment: two incquivalent coupled Hubbard layers (ICHL). Among the new features 
found we mention the appearance of an energy gap in the off-diagonal one-particle Green function. We suggest that 
angle resolved photoemission spectroscopy (ARPES) experiments should be set up to measure this gap foreseen for 
this type of materials (Y 2 Ba 4 Cu 7 Oi5 (247) = YBa 2 Cu 3 7 (l23) + YBa 2 Cu40 8 (l24:)). The results presented in this 
paper, specially the ones of Section |J, represent a first step along these lines, i.e., the use of the sum rules of Nolting 
mTj to determine the physical parameters of the theory. 



II. CALCULATION OF THE PARAMETERS VIA SUM RULES 

The model we study is the Hubbard Hamiltonian Il9|] 



H = KAo c 3<y + —n la n lrf - [ic\ a c lG , (3) 



U 
~2 

where c\ a (ci CT ) are creation (annihilation) operators for particles with spin a. n.i a = cj^c^, U is the local interaction, 
yit the chemical potential (we work in the grand canonical ensemble). We have adopted the Einstein convention for 
repeated indices, i.e., for the N s sites (labelled by i), the z nearest-neighbour (n.n.) sites j and for spin up and down 
(ct = — a = ±1). ti j = —t, for n.n. and zero otherwise. 

In this section we will use the first sum rules or moments of Nolting jlTj applied to the Hubbard model of Eq. (|^) 
to find some conditions on the parameters of the theory. Before applying the sum rules, let us present the spectral 
function, A(k,Lu). It is given by 



A(k,L0) 



5(a) 



1/2 ~ w ) 



sin [<j> - 7r(l/2 - a)] . 
sin0 
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(4) 



if < u> < u! c . In this equation £r = s? — fi e ff, V = u p /u a represents the ratio of the charge and spin velocities and 
0(x) is the usual theta function. The effective chemical potential is defined in Eq. (^). A similar equation can be 
obtain also for the case — lj c < lu < 0. The spectral function (Eq. (||)) has to satisfy the time reversal symmetry, 
condition which leads to 4> = —ira/2, a result already obtained by Yin and Chakravarty |l3| ]. Eq. (||) is assumed to 
be the many-body solution to the Hubbard Hamiltonian (Eq. (ph). 
The first moment (or zeroth order sum rule) M a (k), is given by 



M {k) = 



+ 00 



A(k,u>) dw = 1 



(5) 



Let us say that the condition given by Eq. (g) represents the equal-time anticommutation relation of fermions as 
pointed out by Yin and Chakravarty |13j] . Eq. (|5|) is valid for any fermionic theory, independent of the model used. 
We stress the fact that the fulfilment of Eq. (^) implies that our \PGF describe fermion quasi-particles. Furthermore, 
the area of the distribution is one, i.e., it is normalised. Doing the integration of A(k,u) and using the expression 
given in Eq. (Q), we find 
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(6) 



where B(x, y) is the usual beta Euler function. If we restrict ourselves to regions close to the effective chemical 
potential, i.e., £g « 0, then the normalisation factor, g(a), is independent of k and is given by 



7(«) 



2sin(^; 



(7) 



which reduces to 1 when a — > 0. 
The first order sum rule is given by 



Afi(fc) 



= e% — p + pU = & ; p, e ff = p,-pU 



(8) 



The integral equals to — /i + pf7 is an exact result within the working scheme of the Hubbard model (model 
dependent In Eq. p is the carrier number per lattice site and per spin. We work in the paramagnetic phase, 
namely, p-f = pi = p. Combining Eqs. (||,||), i.e., doing the integral of u A(k,u>) to find the centre of the distribution, 
we get 
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where -F(a, /3; 7; z) is the hypergeometric function. Using an expansion around ^ k as previously, we find that 

2g(a) sin [7ra/2] 



(9) 



?7- (1 -77) Q 



(10) 



which leads to the following relation between the anomalous coefficient a and the spin-charge characteristic ratio 77 

77 — 1 , , 

If we analyse Eq. (|ll|) we can see that as long as0<a<l/2 and < r\ < 1 the only possibility is that a — and 
?7 = 1, which is actually characteristic for the usual non interacting Fermi liquid. As a conclusion, in order to satisfy 
the first order sum rule we have to introduce a new coefficient 7 related to the band renormalization factor. With 
this coefficient the energy £g will be renormalized as 7^. By applying again the first order sum rule, we get 



so the band renormalization factor will be a function of the two previous parameters a and r\. We have to mention 
that the introduction of the band renormalization factor does not affect the time reversal symmetry and the zeroth 
order sum rule. We remark that this relation had not been previously calculated in the literature because only the 
first sum rule (Eq. (|J)) had been used. Now, we are going further and we will apply the second order sum rule or the 
width of the distribution function, A(k,uj), to calculate a. Indeed, what we will find is that this sum rule imposes a 
condition on the three remaining parameters of the theory, namely, w c , a and U. The interaction U is coming from 
the Hubbard model. As the parameter a is not really independent of the Coulomb interaction, our result evidently 
departs from the approach adopted by Tifrea p3] in his Ph. D. thesis and related works [^lj (see the discussion after 
Eq. ©). 

The second order moment, M%{k), is given as 

/+oo 
u; 2 A(k, uj) duj = {e % - p) 2 + 2U{e % - p) + pU 2 = f§ + p{l - P )U 2 . (13) 
-oo 

This is an exact relation and the model dependence is visible thru the local Coulomb interaction. Performing the 
integral of cu 2 A(k, u>) to find the width of the distribution and using the definition of the one-particle spectral function 
given in Eq. (^), we find that in the limit £g « 



2g(a) sin [7ra/2] 



ira a - 



u>i « p(l - P )U 2 (14) 



" Pd-P)f-V • d5) 



a + 2 \0J, 



which gives 



Thus, as we already pointed out, Eq. (|15|) puts a strong condition on the remaining physical variables of the theory. 
We immediately see that a = for U = 0, as it should be. Also, another parameter entering in the constraint is 
the electron number/spin. As < a < 1/2, we see that U max ~ u> c , for p ^ 0; 1. Here we appreciate the difference 
with the work of Tifrea pfj in his Ph. D. thesis where he calculates the superconducting critical temperature, T c , 
in the presence of repulsive local Coulomb interaction. The non-Fermi liquid parameter, a, continues to be a free 
parameter in Tifrea's approach. Before we leave this Section, we say that our energy scales are ordered in the following 
way, T c < A(0) < ujd « u c « D, where A(0) is the superconducting gap at zero temperature, ujd is the Debye 
frequency giving origin to the superconducting critical temperature, T c . 

III. DYNAMICAL AND GLOBAL QUANTITIES 

In this Section we will calculate some dynamical properties of the theory, namely, the momentum distribution 
function, n c (k), in the normal phase at T = 0. Also, we will recalculate the superconducting critical temperature, T c , 
and the superconducting critical interaction strength, A cr , paying duly attention to the presence of the new parameter 
of the theory, 7 = 2/(77 + 1 - 2 ") ( E q- ©)■ 

A. Calculation of n c {st) at T = 

The correlated momentum distribution function is given by the following expression 



exp (u) / 1 ) + 1 



oc 



At T = we have to look carefully to this integral because as long as u> > the exponential function is infinity, 
which implies a zero contribution from the integral. We still have to integrate over the region lo < where the 
exponential function is zero. 

Performing the integrals we end up with 
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(17) 



where f(a, 77) = J Q <izz 1//2 [z + (1 - T])/r)] a 3 / 2 . 



As we see from Eq. (|17|), the correlated momentum distribution function has been calculated close to the effective 
chemical potential. Also, we observe that for a — and r\ = 1 we recover the jump at the chemical potential, as 
it is the case for a Fermi liquid. For a / and rj = 1, this jump has gone away, but the derivative at the effective 
chemical potential is discontinuous. A calculation of the renormalization factor Z, defined as Z — n c (£t) — n c (t;~ ), 
gives Z = 0, a result which implies that our theory is a non- Fermi liquid one. In the other case, a = and r\ ^ 1, we 
obtain the case of a Fermi liquid with spin-charge separation (Z = 1). Another conclusion that we reach by looking 
at Eq. ([It]) is that the parameter 7 = 2/ (77 +1 — 2a) modifies the results of Yin and Chakravarty [|l3| in the following 
way: in order to study non-Fermi liquid systems we have to consider that the frequency cutoff is effectively smaller, 



i.e. 



uj c — > (1 - a)uj c 



B. Calculation of T c and A c in an s— wave superconductor 

In the following we will restrict ourself to the study of the non-Fermi liquid system (a =/= 0, r\ = 1) where the Green's 
function according to our previous calculations is given by 

G (k,u) = 7 ' , — (18) 

with 7 = 1/(1 — a). The order parameter equation in the framework of the Gorkov equations is given by: 

1 = GoHkMGo\-k,-iu; n ) - |A fc H 2 (19) 

where (3 — 1/T. The critical temperature will be obtained from Eq. (|l9| ) with the condition — > 0, and the difficult 
problem of evaluating the sum over the Matsubara frequency will be solve by using a contour integral similar with 
the one used in Ref. |20|]2l|]. In the limit j3u>n <C 1 the critical temperature can be obtained exactly as 



rp2a _ 1 



C(a) 



D{a) (^uj D ) 2a 



g 2 (a) XA(a) 



(20) 



where A = 1/N(0)V and A(a), C(a) and D(a) has the same meaning as in the paper of Muthukumar et al. [El]. 
We have to mention that our critical temperature is different from the one obtained in Ref. pl[ , and include the 
renormalization factor g{a) and the band renormalization factor. As a result we can see a decrease of the critical 
temperature due to the effective Debye frequency — jojd ■ We also obtained a modified critical coupling constant 

Xcr = ?fe) (^) A(a)D(a) (21) 

which by the same reason of an enhanced effective Debye frequency seems to be enhanced. In Fig.l we present a plot 
of the critical cou plin g constant versus the non-Fermi parameter a. As we can see in the limit a — > we recover the 
usual BCS result @ (A cr = 0). 



IV. THE ISOTOPE EFFECT FOR NON-FERMI LIQUIDS 



The isotope effect exponent is an important physical parameter since, in the low temperature superconductors 
(LTSC), it was used to determine the origin of the pairing mechanism. In the HTSC the isotope effect has been 
widely discussed, in particular to detect the influence of the phonon degrees of freedom on T c . As Kishore [^5| points 
out, the isotope effect in the cuprates is sensible to various factors like: 1) the form of the free density of states, 2) the 



Coulomb interaction, 3) the carrier concentration, 4) presence of impurities, 5) anharmonicity and 6) the symmetry 
of the order parameter, among others. From our Eq. (|l5|) we can account for the first three dependences. Here we 
will calculate T c and the isotope exponent, a', for a s-wave superconductor. Thus, our expression for a' could be 
approximately applicable to the cuprates to fit the data. We believe that the numerical value of global quantities, 
like a', do not depend too much on the symmetry of the order parameter. It appears that La2-xSrCuO^ presents a 
s-wave symmetry order parameter. So, it is the natural candidate to apply the ideas worked out here. On the other 
side, BiiSriC aCuiO% [|26| holds a d-wave symmetry order parameter. 
The isotope effect exponent is defined as T c cx M~ a . Then, it is given by 



<91nT c 
"dlnM 



(22) 



where M is the isotope mass. Another relation that we will use is the relation between the Debye frequency and the 
isotope mass, namely, u>d oc A/ -1 / 2 . Now, Eq. (EOT) can be re- written as 



~ = g 2 (a)A(a)F(cj D ,a,T c ) 
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(23) 



Deriving Eq. (E3J) with respect to M we get that 



( 9F 



2T r _ 



(24) 



Let us mention that Eq. (J2J) is valid for the case that a ^ /(cj£>). The case a — /(wo) will be discussed later on 
(Eq. (p6|)). Performing the partial derivatives, we come down to the following expression 



1 - 
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(25) 



From Eq. (pjJ) we gain the BCS case, i.e., a 1 = 1/2 when a = 0. In Fig. 2 we plot the isotope coefficient versus the 
non- Fermi parameter a for different values of the coupling constant. As we expected there is deviation for the usual 
BCS result as a increases. Also we mention that the isotope coefficient depends on the coupling constant, a result 
similar with the one obtained for the 2A(0)/T C ratio |p0|| . 

While leaving this Section, let us mention that if a would depend on lod, then the isotope effect exponent should 
be obtained from the following expression 



UJ D 
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(26) 



We could make a = fiuio) by choosing ui c oc ujjj, as Tifrea has done in his Ph. D. thesis Q. This contribution is 
important to get full agreement with the experimental data, namely, that a' ~ at optimal doping. We argue that 
the lowering of a' from the BCS result should come from the second contribution in the numerator in Eq. ( |2q ) . We 
leave this task for the future. However, the result found in the present section should be applicable in the underdoped 
regime of the HTSC's. 



V. CONCLUSIONS AND OUTLOOK 

We have applied the first three sum rules of Nolting Jl?]] for the NFLGF of Eq. (Q), which is assumed to be a 
solution of the Hubbard model (Eq. (||)) for frequencies |a>| < uj c . This NFLGF is anomalous due to the presence 
of the non-Fermi liquid exponent, a. Due to the requirement that the first order sum rule of the spectral function 
be satisfied, a new parameter has to be called for. 7, the so called mass renormalization factor plays an important 
role in the theory. Due to its presence we have recalculated some results found in the literature and pointed out the 



role of 7. Also, a new quantity, the isotope effect exponent or a', has been calculated for the first time for non-Fermi 
liquid systems, a' reduces to the BCS case when the Coulomb interaction is zero, namely, for a = 0. The fact that 
a oc U 2 comes out from the application of the second order sum rule to the spectral function. Due to Eq. (|l5|), a 
is a model dependent parameter, a result found previously in the literature by Wen ||. According to our view, T c 
decreases with a p0| ] and a (x U 2 , then we can conclude that the local Coulomb repulsive interaction is detrimental 
to superconductivity. This result has been found in other approaches p8p9| ]. 

We would like to point out that our spectral function (Eq. (Q)) does not reduce to the well known Dirac delta 
function when a = 0, 77 = 1. In order to re-obtain the usual delta form of the spectral function we have to set a = 
and 77 — 1 in our starting Green function. Even without doing that for a = and 7/ = 1 the spectral function (Eq. ||) 
will satisfy the usual scaling relation for a Fermi liquid system A(A[k — fop], Aw) = A _1 yl([fc — fcp],u;). We used this 
property in some parts of the paper and we have recovered known features, like: 1) the jump at [i of the momentum 
distribution function, n c (k), which is a Fermi liquid behaviour; 2) the superconducting critical interaction is zero, 
namely, X cr = 0, and 3) the BCS isotope exponent, a' — 1/2. 

In Appendix |A| we have applied the formalism developed in this paper to two inequivalent coupled Hubbard layers 
(ICHL) extending the calculation of Yin and Chakravarty jl3| for equivalent planes. Hildebrand et al p7|] have 
applied the FLEX formalism to ICHUs. In particular, our expressions for G\\{k, ito n ), G 2 2{k 1 iu n ) and G\ 2 {k, iu) n ) — 
G2i(k,iu> n ) are different from the ones of Ref. jl3| due to the presence of a shift in the effective chemical potential, 
i.e., fJ.i.eff = M2,e// + S. Experiments should be designed to detect the results found in our work (See Appendix [A|). 
In particular, the presence of the theoretical gap in off-diagonal one-particle spectral function, A\2{k,ui), calls for 
ARPES experiments to be performed in these materials |3(| . 

We have assumed that the non-Fermi liquid parameter, a, is independent of k. This seems not to be the case 
as Meden |3^] has pointed out. He concludes that the asymptotic behaviour of the one-particle Green functions of 
Luttingcr liquids at large space time distances is not universal. Namely, along certain directions the exponent of the 
asymptotic power law is not given by the Luttinger liquid parameters. Due to this consideration, a could depend on 
k. This possibility is outside the scope of the present paper. Among one of the possibilities we would like to explore is 
the calculation of pressure effects |3^| in some HTSC materials, when correlation is important, i.e., a ^ 0. Another 
possibility worth considering, for the case of inequivalent coupled Hubbard layers, is ot\ ^ aa, as the two non-Fermi 
liquid parameters on the two planes. This approach is likely more demanding than the one followed in the paper. 



We leave for the future the self-consistent calculation of p\ and P2 on the two layers (See Eq. (A5)). On the other 
hand, the presence of the mass renormalization parameter, 7, will modify the zero temperature order parameter, 
A(0). However, its calculation has been left out of the present work. Of course, with A(0), we could calculate the 
ratio 2A(0)/T C = f(U) 7^ 3.5, seeing its dependence on U (or a). Another aspect which should be addressed in the 
future is the superconducting properties of two inequivalent coupled Hubbard layers. In this case, our mean field 
Hamiltonian becomes a 4 x 4 matrix in the Nambu formalism, giving rise to pairing on the planes and perpendicular 
to them. We would like to end by saying that the nature of the superconducting transition is strongly related to how 
anomalous (non-Fermi liquid like) the normal state spectral function is, and as such, is dependent upon the doping 
level ]33| ] . The anomalous properties of the normal state spectral function are visible in the underdoped regime of the 
HTSC. 
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APPENDIX A: NON-FERMI LIQUID ONE-PARTICLE GREEN FUNCTIONS FOR TWO 
INEQUIVALENT COUPLED HUBBARD LAYERS 

In Ref. p7| , Hildebrand et al have studied the case of two inequivalent coupled Hubbard layers for the case of 
Y 2 Ba 4 Cu 7 15 {247) = YBa 2 Cu 3 7 (123) + Y Ba 2 Cu i O ii {l2A). We will re-study this system using NFLGF for each 
one of the layers with the purpose of finding new theoretical consequences of this assumption. For example, we are 
going to compare with the results of Refs. ]l3| , p7|j . The Hamiltonian of the system is given by 

H = ( G^ l {e^,iu n ) t ± \ 



where t± is the coupling matrix element between the two inequivalent Hubbard layers and Gj(e^,iuj n ), j = 1,2, are 
the normal state non-Fermi liquid one-particle Green functions of the layers. They are given by 



9(a) 



1 



^ctg±i7ra /2 



(lOJ n - 7^. gj 



(Al) 



with rj 1 g = eg — fJ- e ff an d ^2 fc = 7 ?i £ ~ ^' ^' e '' we nave included a shift between the two effective chemical potentials. 
A simple calculation allows us to calculate the diagonal and off-diagonal one-particle Green functions. They are 



Gjj(e^,iu n ) 
Gi 2 (ej:,iuj n ) 



G- 1 (e^,iuj n ) 



(G 1 1 {e^iu n )G 2 1 (s£,iuj n ) - t 2 ± ) 
tx 



(G 1 1 (e^iu n )G 2 1 (eg,ico n ) - t 



; j = l,2 ; 1 = 2; 2 = 1 
, G 2 i(e^iuJn) = Gi 2 (ej,iw„) 



(A2) 



The excitation spectrum is determinated by the roots of the denominator of Eq. (|A^) . There are branch points at 
lu = 777 . g. Let us assume that 77 1 g > r\ 2 g > 0. Thus, for > 0, we must divide the complex plane in six regions 
named I-II-III-IV-V-VI, as shown in Fig. 3. a. The poles are given by the solutions of 



e i™ul a g- 2 (a) (u -77^ g + i0+) " L -(3r] 2 j : + i0 + ) " =t\ e 2i " * ; n an integer , 



(A3) 



Following the analysis performed by Yin and Chakravarty |13| we conclude that, for r] 1 g > r\ 2 g > 0, there are solutions 
only in the regions denoted by I and V. In the case that 0<77 1 g> 77 2 ^<0, still for oj > 0, we divide the complex 
plane in four regions I— II— III— IV as shown in Fig. 3.b. Similarly to the previous analysis, we conclude that there 
is solution in region / only. The analysis is similar for the case rj 1 g < rj 2 g < 0. 
The poles are localized at 



7(%E + V 2 .k) ± 



^1.2 



V 2 1 



■At 2 , 



,eff 



1/2 



ff(a)tj 



(A4) 



From Eq. ( A4) we recover the case of equivalent coupled planes of Yin and Chakravarty J13| by making r\ x g = r\ 2 g. 

Now, we are in a position of calculating the one-particle spectral functions, namely, Ajj(fc,oj), with i,j — 1,2. The 
results are (with jr] 1 g > 77^ g) 



An(fc,w) = 



where Xa 



A 2 2(fe,Cj) = 



(l + (X 2 Xi)''°) 

sm(7ra/2) A 2 ^ Q J (i+(*2 x 1 ) 2 < 1 -°'-2co S ( TO )(x 2 
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if w < 7^ < rmi,i or 7»7 2l S < 7%,£ < ^ 

if 77y 2 g < oj < 7 7 7l g 

if w < 777^ < 77 ?1 g or 7 77 2 g < 777^ < w; 

if 777^ < w < 7 77 a 



A 12 (fc,w) = A 2 i(fc,w) 



sin(7ra) (X\ X 2 ) 
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if 7771 < cj < 7772 



In Figures 4 we show the behaviour of the diagonal and off-diagonal one-particle spectral functions. From Figure 4. a 
we observe that the symmetry of the diagonal spectral function, An (k, oj), is lost around the frequency u> = r yq 2 g = 0.9. 



To realize the new feature, we refer the reader to the respective figure (i r n(x)) in Ref. Jl3|]. Also, the symmetry of 
^22(2^) in Ref. |l3f| is lost around to = 777. g = 1.2. In consequence, ARPES experiments shining light on layer 1 
or layer 2 should detect these fine details. Namely, we will have two symmetry breaking, 1 — ► 1 and 2 — > 2. Now, 
from Figure 4.b we see that there is an interval of energy where A 12 (fc,a->) = 0. This feature was not obtained for 
the off-diagonal one-particle spectral function for equivalent planes (See F\2{x) in Ref. [^3[ ). Another aspect which 
we would like to point out is the fact that in the one-particle spectral functions there is not an unique variable to 
describe the data. We have two relevant energy scales in the problem, for each value of k. For the type of materials 
we are studying it seems natural to find, on theoretical grounds, a gap in the energy spectrum. ARPES experiments 
should be designed to detect the gap found from our expressions [|(j| . We mention that the gap found in A12 (k, oj) 
depends on the relative carrier number in the two inequivalent coupled Hubbard layers. Therefore, the presence of 
a shift between the two effective chemical potentials produces new theoretical results for the spectral densities. In 



particular, non-Fermi liquid quasi-particles aquire a more complex structure (Eq. (A4)) than the case of equivalent 
planes. 

In order to calculate the carrier number per site per spin per plane, at T — 0, we have to perform the following 
integrals 



1 



+ D r 



-D 



Pi(2) = 7 / ds I N(e) A 11(22 )(£,w) du , (A5) 



which have to be computed numerically. Work is in progress to solve Eq. ( |A5| ) self-consistently. N(e) is the 
uncorrelated density of states. To conclude this Appendix, we say that our results generalise the case of Yin and 
Chakravarty |lj|]. With respect to the results of Ref. [^Tj, they have not discussed the two one-particle spectral 
functions, An(k,u>) and A 12 (fc, oj), around the effective chemical potential. Most likely the results of Ref. p7| are 
valid for large values of frequencies. 
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FIG. 2. The isotope coefficient versus the non-Fermi parameter a for different values of the coupling constant ((1)-A = 0.33, 
(2)-A = 0.66. (3)-A = 1). 
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FIG. 3. (a) For uo > 0, ~/r] 1 g > 7?? 2 g and 77^ 2 jg > 0, we have divided the complex plane in six regions, showing the branch 
points (jrj 2 g and jrj 2 g) and branch cuts. The branch cuts are given by the dash lines. The solutions, for w > 0, are in regions: 
1) /, namely, n = 0; and 2) V. (b) For a; > 0, 7^ g > 77) 2 g, and jj 1 g > 0, ij 2 g < 0, we have divided the complex plane in four 
regions, showing the branch points and branch cuts. For ui > the solutions are in region I. 
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FIG. 4. (a) The diagonal one-particle spectral function, Au(k,u>), vs u), for some values of r YH 1 r and 7?7 2 r, namely, 
jr) l g = 1.2 and 7?7 2 g = 0.9. Compare with Fu(x) of Ref.[13] We notice that the symmetry around to = 0.9 is lost. This is a 

realization of inequivalent coupled Hubbard layers. In similar form, A 2 2(k, u), vs u should not be symmetric around uj = 1.2, 
according with the given parameters. (b)The off-diagonal one-particle spectral function, A 1 2(k,uj), vs u). Same values as in 
Figure 4. a. Compare with F\2{x) of Ref.[13]. Now we have a gap in the energy spectrum betweem lj — 0.9 and to — 1.2. This 
is also a consequence of having inequivalent coupled Hubbard layers. 



